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I. INTRODUCTION 
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(Dated: October 3, 2008) 

Domain walls in a tetragonal chiral p-wave superconductors with broken time reversal symmetry 
are analyzed in the framework of the Ginsburg-Landau theory. The energy and the jump of the 
magnetic induction on the wall were determined for different types of walls as functions of the 
. . . parameters of the Ginzburg-Landau theory and orientation of the domain wall with respect to the 

crystallographic axes. We discuss implications of the analysis for Sr2Ru04, where no stray magnetic 
I fields from domain walls were detected experimentally. 
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There were a number experimental evidences of the broken time-reversal symmetry (TRS) in the unconventional 
superconductor Sr2Ru04 [l], 0, Q. It was suggested that this phenomenon was connected with the p-wave Cooper 
Q ' pairing with the wave function in the momentum space proportional to px + ipy (chiral p-wave superconductivity) 
In chiral p-wave superconductors the spontaneous magnetic flux must be present near domain walls or sample 
5^ ■ boundaries. However, Kirtley et al. Q have not detected any stray magnetic fields, which the spontaneous magnetic 
CIh' flux should produce above the sample surface. This put in question the scenario of the chiral p-wave pairing and 
^ . stimulated theoretical investigations of the problem. In particular, the relation between the stray fields and the 
^ ' magnetic flux, which appears near the domain wall (DW) in the bulk, was derived 0, Q • Another challenge for the 
, theory was to find the spontaneous magnetic flux itself. It is determined by the product of the London penetration 
' depth and the jump of the magnetic induction on the domain wall. The latter is on the order of the first critical 
, , magnetic field Hci and was already analyzed in the past for particular types of the DWs in chiral superconductors 

^ \ There are two possible explanations why the experiment could not detect the stray fields above the sample surface. 
O , The first one is that there is a domain structure with a period so small that stray fields decay very fast in space 
at a distance small compared to the distance of the experimental probe from the sample surface. There were some 
experimental evidences of domain structure in Sr2Ru04 [l^, [11], though Xia et al. [3] did not reveal any domain 
structure at studying the Kerr effect. The theory predicts that in superconductors with broken TRS usual ferromag- 
^ , netic domains, which depend on the size and the shape of the sample (extrinsic domains), cannot appear p^ . but 
another type of domains, which decrease the bulk magnetostatic energy at the expense of destroying the Meissner 
state, becomes possible [H, [l^ . The size of these domains is roughly of the order of the London penetration depth A 
r. and does not depend on either shape or size of the sample (intrinsic domains, see discussion in Ref. ,17,). The analysis 
- , • ■ of intrinsic domain structure in superconductors with broken TRS was recently extended on the case of finite external 
magnetic fields, and the equilibrium magnetezation curves in the state with intrinsic domains (cryptoferromagnetic 
state) were found theoretically [isj . The intrinsic domain structure may lead to serious suppression of stray fields 
, above the sample, though definite quantitative conclusions are not yet possible because of the absence of essential 
information on Sr2Ru04 (magnetic crystal anisotropy, e.g.). Therefore, the other possible explanation of the negative 
result of Kirtley et al. Q must be also considered: a very small spontaneous magnetic flux penetrating along DWs in 
Sr2Ru04. This problem is investigated in the present work. 
?H ' The magnetic flux near DWs originates from the intrinsic orbital moment of Cooper pairs (orbital ferromagnetism) . 
As was demonstrated in Ref. [l^, in the case of orbital ferromagnetism one should not use the Landau-Lifshitz theory 
[20| , since in this case one cannot define local spontaneous magnetization. Therefore one should rely on magnetization 
currents generated by the orbital moment of Cooper pairs, which cannot be reduced to the curl of any vector moment. 
The total magnetization current along the DW leads to the jump of the magnetic induction on the DW, which 
eventually determines the magnetic flux connected with the DW. In the present work we extend previous theoretical 
investigations of the magnetic flux around DWs in tetragonal chiral p-wave superconductors analyzing conditions 
for appearance of different types of DW structure depending on the strength of crystal in-plane anisotropy and DW 
orientation with respect to crystallographic axes. The analysis was done using the Ginzburg-Landay (GL) theory, 
and the magnetic induction near DWs depends on the parameters of the GL theory. For the parameters obtained in 
the weak-coupling limit the magnetic induction near DWs is scaled by the first critical magnetic field Hd differing 
from the latter by a numerical factor. But for other values of the GL parameters the magnetic induction near DWs 
can be much smaller and even become negative with respect to the intrinsic magnetic moment of Cooper pairs. This 
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stresses again that the orbital moment of the Cooper pair in the chiral p-wave state does not lead to any definite local 
magnetization. 

The practical outcome of the presented analysis is that one can explain very weak stray fields from DWs if the 
GL parameters differ from those calculated in the weak-coupling limit. The final conclusion on the reason for the 
absence of detectable stray fields around Sr2Ru04 samples is possible only after these parameters are found from the 
experiment or from the strong-coupling theory. 



II. MODEL 



The unconventional superconductor Sr2Ru04 belongs to the tetragonal crystallographic symmetry group -D4/1. 
Considering the p-wave state for this material they usually believe that strong crystal anisotropy keeps both spin 
and orbital momentum of the Cooper pair parallel to the c axis. Then the p-wave state corresponds to the two- 
dimensional representation = {kxzjtyz}, and the order parameter {ijxkx + r]yky)z is fully described by a complex 
two-component vector rj = iVx,Vy) K 5|. Then the GL free energy density is 



+K3[iDxVxT{DyTjy) + {DxVx){Dyr,yy] + K4{DxVyT{DyTjx) + {DxVy){DyT]xT] + K^{\D,T^x? + \D,Tiy\^). (1) 

Here Di — di — i{2e/ch)Ai is a covariant derivative and B = V x A is the magnetic induction. In the case of full axial 
symmetry in the plane (3^ — and Ki = K2+K3+K4, In the BCS theory (weak-coupling limit) Ki/3 = K2 = K3 = K4 

The complex components of the two-component order parameter r/ can be represented in a form 

Vx = Vo cose exp(i0), 
rjy^rjo sine exp[i{<j) + Q], (2) 

where angles (phases) 6, (p, and ( are coordinate dependent functions in general, and t]q = \ri\. In the focus of 
the present study is the chiral state with the two-fold degenerate ground state 77 = 2~^/^77o exp (i0)(I, ±i), which 
corresponds to e = 7r/4 and ( = ±tt/2 and is reahzed if /32 > and /3i > P2 — Ps/^: > 0. The order-parameter 
amplitude, 

^0 = 4(/?i-^V+/33' 

is determined by minimization of the GL free energy. The chiral state has a nonzero orbital moment I — ir] x rj* = l^z 
with the only component 

Iz = iiVxVy - VyV*x) = sin2esinC. (4) 

Other components of the orbital moment are suppressed by strong crystal anisotropy. Though they are necessary for 
orbital magnons , they can be ignored in our analysis of DWs. 

Substituting the general expression for the two-component order parameter, Eq. into the free energy density J- 
[Eq. ([T])], one can derive the expression for the superconducting current: 

e 

TO 

which consists of two contributions. The transport current j'^*''-' — env is determined by the gauge-invariant superfluid 
velocity 

V - T5^V(/) - — A, (6) 
Zm mc 

where n is the superfluid electron-density matrix with the components 

rixx = ^VoiKi cos^ e + K2 sin^ e), 

8m 2/ 2 2 \ 

Uyy = ^riQ(K I sm e + ^2 cos e), 

nxy = riyx = -^Vo {K3 + K4) sin 20 cos C, 



j = ^a^^ = j(f)+j(™), (5) 



-j^VoK5- (7) 
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The second contribution to the current is the magnetization current with the components 

= ^r]l[2K2 sin^ 6 d^C + 2{Kz cos^ 8 + ^4:4 sin^ 8) sinC dyQ + K3 sin2e cosC dy(], 
jj") = yTygpifi sin^ e dyC + 2{K3 sin^ 8 + i^4 cos^ 8) sinC 9:^6 + i^4 sin 28 cosC d^C], 

jt"^ = ^77o^5sin2 8a.C. (8) 

In the following we shall neglect the terms oc K5 considering only DWs parallel to the axis z when parameters do not 
vary along this axis. We also assume that the London penetration depth A is much larger than the thickness of the 
DW, which is on the order of the coherence length ^. This allows to ignore Meissner currents studying currents inside 
the DW. 

Inside domains 8 = 7r/4, and n^x — nyy — (4m/fi^)?7g(i4'i + K2). Then the London penetration depth and the first 
critical magnetic field are 

The magnetization current cannot be presented as j*^™-' — (l/c)V x M, which is the main assumption of the 
Landau-Lifshitz theory. Therefore for chiral superconductors the Landau-Lifshitz theory is not valid, and the local 
magnetization M cannot be defined [19]. Whereas in the Landau-Lifshitz theory the jump of the tangential component 
of the magnetic induction on the DW is given by the universal relation 5B = 8ttM independently of the DW structure, 
in chiral superconductors 6B does depend on the DW type and on the DW orientation relative to crystallographic 
axes. This will be demonstrated in the next section. 



III. DOMAIN WALLS 



The two degenerate ground states correspond to two possible directions of the orbital moment Iz — parallel or 
anti-parallel to the c axis (the axis z). This leads to the existence of domains with Iz — ±1 separated by domain 
walls. Two types of DWs are known. The DW of the type I is characterized by a gradual change of the phase ( = C(r) 
at constant 8 — 7r/4, i.e., the absolute values of the components \rix\ = {riyl remain constant inside the DW We 
address this type of the DW as C— wall. In the DW of the type II the relative phase ( — —tt/2 remains fixed, and the 
transition from 1^ = —1 to 1^ = +1 is realized via variation of the angle 8 so that in the center of the DW one of the 
components rj^ or ijy vanishes. This type of DW was considered by Sigrist et al. [lo| assuming that the component of 
ry, which does not vanish, remains constant. In the present work we assume that inside the DW the order parameter 
rj rotates in space at fixed order parameter modulus \ri\ = rjo. This assumption is justified if /3i ^ p2,l33- Further we 
address this type of the DW as 8— wall. However, the difference between the two types of DW depends on a choice of 
the coordinate system: rotating the in-plane coordinate frame xy, in which the components of 77 are defined, through 
the angle n/A (45°), the ^-wall becomes a 8- wall and vice versa. 



A. Variation of relative phase: (^-wall 

Let us consider a DW, which is characterized by continuous variation of the relative phase C = C(r) at fixed angle 
8 = 7r/4. The transformation of the Iz = —1 domain to the Iz = +1 domain is possible via counterclockwise rotation 
of C from ~tt/2 to tt/2 or clockwise rotation of C from ~tt/2 to —3tt/2. The structure of this DW, its surface energy, 
and the jump of magnetic induction depend on its orientation with respect to chosen crystallographic axes (xy). 
Assuming the arbitrary orientation of the DW plane in the xy-j>\a,ne we transform the gradient terms in Eq. ([T]) to 
the new coordinate system (x'y') connected to the DW: the axis x' is normal and the axis y' is parallel to the DW: 

dx — cos il>dx' — sin ipdy> , 

dy — sintljdx' + cosi/jdy'. (10) 

The transformation is a counterclockwise rotation of the original coordinate axes xy around 2;-axis through an angle ip 
about the origin. Note that only the gradients were rotated, the order parameter 77 being untouched and characterized 
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by the components rj^^y in the original coordinate system xy. After the transformation the free energy density with 
respect to the energy of the uniform chiral phase is 

AT= ^^(1 +COS2C) + ^^^{{cosi:D^,(j) - smi;Dy,(j))^ + [smtP{D^, cj) + d^>C) + cos ij{Dy,(t) + dy>C)?} 



{{s'mtpD^xj) + cosipDy>(l)f + [cos ip {D cf) + d^'C) - smipdy'{(j) + C)]'^} 



+7^3770 cosC(cos — sin^Dyi (j))[sin'ip{Dx' (j) + dx'C) + cos + dyiQ)] 

+Ki'ql cos({sin^pDx'4i + cosipdyi(j))[cosip{Dj:'(p + dx'C) - sin-tp^DyKj) + dy'C). (11) 



The relevant components of the superfluid electron-density matrix are 



Am 



Tlx'y' — '^y' x' 



ril [Ki +K2 + [Ks + if 4) sin 2iP cos C] , 



4m 



Vo {^3 + ^i) COS 2ip cos C- 



(12) 



The magnetization currents are defined by 



.(m) 



-772 [Ki + (Ki - K2) cos 2V + (Ks + K^) sin 2^/. cos C]9x'C, 



4"^ = {[-^3 + Ki + {K^ + Ki) COS27A] cosC + (i^i - K2) sin2V'}a,<. 

,(m) 



(13) 



The current normal to the DW plane must vanish: jx' = enx'x'Vx' + jx' —0- This gives the following expression for 
x'-component of the velocity: 



2m 



h 

4m 



(Ki- K2)cos2il; 



K1+K2 + (if 3 + K4) sin 2ij cos C 



dx'C- 



After its exclusion the total current jy' parallel to the DW is 

{Ki + if 2) sin 2iP + {K3 + is:4) cos C 



-(if3-i^4)cOsC+(ifl-if2 



dx'C. 



K1+K2 + {K3 + Ki) sin 2ip cos C 
Using Eq. psp and the Maxwell equation one can calculate the jump of the magnetic induction on the C-wall: 

'tt 2cos^2?/; 



87re 



5i3c(V) = -(4vr/c) / dx' jy = — ry^ \ {K3 - if4) T 



K,-K2 
sin 2-0 



2 (l-q2)i/2 



arctan 



1T9 



l±g 



where 



g3 + j^4 . 

'^=ifrT^"^'^- 

The upper and the lower signs correspond to C rotations from — 7r/2 to tt/2 or to —3tt/2 respectively. 
After the exclusion of the transport velocity the free energy density becomes 

AT=^il- COS27) + (ifi + K2)v'jil)^, 



(14) 



(15) 



(16) 



(17) 



(18) 



where 



/(t) = 1 T 9sin7 - 



Ki- K2Y cos^ 2iP 



if 1 + if 2 / 1 ± g sin 7 



(19) 



Here we introduced the angle 7 equal to7 = C + ''''/2or7 = — (^"""/S for C-rotations from — 7r/2 to tt/2 or to — 37r/2 
respectively. So across the DW 7 varies from to tt. 

The structure of the (^-wall is determined by the Euler-Lagrange equation obtained by variation of the free energy 
with the density (fTS]) with respect to 7: 



f{i)di'i+\f{indx'i)'-'^ 



0, 



(20) 



5 




FIG. 1: Shown are the angular dependence (a) of the jump of reduced magnetic induction (5_Bj x 3hc/32TverioKi and (b) the 
reduced surface energy cr^ x 3/47;q(2/32^i)^^'^ of the DW for two configurations, when C, increases from —n/2 to 7r/2 (curve 1), 
or decreases from — 7r/2 to ~3n/2 (curve 2) with increase of x' . 



where the scale 



1 fKi + K2^^'^ 



^-2S(^I m 



determines the thickness of the DW. The first integration of Eq. (|20p yields 

L — cos 2 



n9 1 ^ COS 27 



After this one can find the surface energy of the DW: 

<^d^)=vl{^^^^^ [ dl^^l^). (23) 

Figure [T] shows the plots of the magnetic induction jump J_B^ and the surface energy as functions of the angle 
between the DW and one of the crystallographic axes for the weak-coupling case when K2 = = K4 = Ki/3. 
The plots for the clockwise rotation {( varies within DW from —tt/2 to — 37r/2, curves 2) are obtained from the 
counterclockwise rotation (C varies within DW from — 7r/2 to 7r/2, curves 1) by reflection ip —tjj. The angular 
dependences have two extrema at ^ — ±7r/4. The spatial distribution of the magnetization currents along the DW 
for different orientations of the DW is shown in Fig. [2] 

Let us consider some important particular cases. The case ip =0 corresponds to the DW parallel to one from chosen 
crystallographic axes. The surface energy and the jump of magnetic induction in this case are 



ac(0) = 2%^ 



3 12P2K1K2 



K1+K2 
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X'/A 

FIG. 2: (color online) The coordinate dependence of the reduced superconducting current j^i x 3^''^ft/2"^''^er;o(-R'i/32)^''^ at 
different orientations of the DW (|<^| < 7r/2). 



This agrees with the results of Volovik and Gor'kov [9|, who assumed that = K4. In the case ^ = there is 
degeneracy between DWs corresponding to two senses of rotation of the phase the surface energies and the field 
jumps for two cases coincide. 

Another particular case is the angle V' — 7i'/4, where the plots in Fig. [1] have extrema. In this case according to 
Eq. (jl6p the jump of magnetic induction is 



SB, 



87re 



(25) 



and is different for both configurations. The DW surface energy can be obtained from Eqs. (fT9|) and ((23|) : 



TT 

4 
3F2 



= Vo 



P2{Ki + K2 



■ / ^7 sin7\/l =F 9sin7) 
Jo 



11 13 



nq 



1 3 



where 



q ■ 



K1+K2' 

and pFg{ai, ...,ap;bi, ...,bq; z) is the generalized hypergeometric function [22 



(26) 



(27) 



B. Rotation of order parameter: 0-wall 



Let us consider now the DWs with rotation of the order parameter rj in the configuration space. Inside this domain 
wall the phase 8 rotates from — 7r/4 to +7r/4 (counterclockwise rotation) or from — 7r/4 to — 37r/4 (clockwise rotation). 
We address this type of the DW as 0-wall. The whole analysis, which was performed above for the C-wall, can be 
repeated for the 0-wall. This yields the expressions obtained from those for the ^-wall by the following substitution: 
C 28, ip ^ ip + 7r/4, Ki ~ K2 ^ + K4 and (32 ^ 02 — 4/33. In particular, the final expressions for the 
magnetic-induction jump and the DW surface energy are: 



SBeW 
where 



-(47r/c) / dx'jy 



c { cos 2^p 



2-^^"''^ (l-g2)i/2 



arctan 



l±q_ 



K1-K2 _ 



(28) 



(29) 
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and 



where 



ae{iP) = ?7o [2/32 (i^i +X2)] / d^sm-f^/W) 



fil) = 1 T 



K1-K2 
K1+K2 



cos 2'ip sin 7 



sin^ 2-0 



i^i + K2 J 1 ± (? sin 7 



(30) 



(31) 



Here 7 = 28 + 7r/2 or 7 = —20 — tt/2 for 8-rotations from —it/ A to 7r/4 or to — 37r/4 respectively. The DW surface 
energy has extrema at -0 = where the magnetic-induction jump is 



5Be{Q) = —r,l 



(32) 



IV. DISCUSSION 



Let us discuss consequences of the obtained analytical and numerical results. We start from the case of full axial 
symmetry in the xy plane, when Ki = K2 + K3 + K4 and P3 — 0. In this case the -0 dependencies for the and the 
0-walls are identical except for the shift of tp by 7r/4. Naturally the choice of crystallographic axes does not matter in 
this case, only their orientation with respect to the DW being important. The minimum of the surface energy for the 
(■-wall corresponds (see Fig. [T|d) to the angle tp = ±7r/4 (45°) between the DW and the axes for 77. However, choosing 
for 77 the coordinate system connected with the DW, the DW becomes a 9-wall. 

Let us consider now the tetragonal symmetry with 0. Without any loss of generality one may assume that 

< 0. Indeed, if /Ja > with 7r/4-rotation of the crystal axes one obtain the 4th order terms with P2 — Ps/A: instead 
of /32 and — /Ja instead of P3. Thus after this rotation the new is negative. At our choice of crystallographic axes 
the (-wall has a smaller surface energy with the minimum at ip = ±7r/4. Thus the stable DW is a (-wall with respect 
to the order parameter axes, but is a 0-wall, if one use the order parameter rj defined in the coordinates related to 
the stable DW. The stable configuration of the DW corresponds to the maximum jump of the magnetic induction 
given by Eq. (j25p with the lower sign before the second term. On the other hand, the jump is minimal and is given 
by Eq. (|24p if the DW is parallel to one of the two crystallographic axes. It is interesting to compare these extremal 
values with the first critical magnetic field given by Eq. ^ : 

SB„,a^, ^ tt{Ki ~ K2) SB,mn ^ AjKiKj - K2K3) 

Hci (ifi -f if2) ln(A/0 ' Hci {Ki+K2fH\/iy ^ ' 

For the values A = 190 nm and ^ = 66 nm ^] and assuming the weak-coupling relations Ki/Z — K2 — K3 — K4, 
these ratios are 1.48 and 0.48 respectively. 

Since now there is no freedom for the choice of the axes for the order parameter 77 there is a force, which tends to 
orient the DW at 45° to the crystallographic axes. This force may compete with other forces on the DW, e.g. , those 
from the sample shape. This would result in various configurations of DW not necessary those dictated by the bulk 
tetragonal anisotropy. 

Whereas in the Landau-Lifshitz theory of ferromagnetism the magnetic-induction jump is given by the universal 
value 87rM, where M is the local spontaneous magnetization, in our case this jump depends on the structure and 
the orientation of the DW. Moreover, if one tried to introduce formally the effective magnetization via the relation 
M = SB I Stt the latter has no straightforward connection with the orbital moment 1 of the Cooper pair and even can 
have a negative sign with respect to the intrinsic magnetic moment of the Cooper pair. 

Returning back to the question why they could not detect stray fields generated by the magnetic-induction jumps 
outside Sr2Ru04 samples @ our analysis demonstrates that at the present stage the theory cannot make definite 
predictions on the strength of these fields without reliable information on the parameters in the GL theory. In 
principle, one could suggest that the DWs in the experiments were not in the ground state, or the GL parameters are 
essentially different from their weak-coupling-limit values. But further experimental and theoretical work is needed 
to check these suggestions. 



V. CONCLUSIONS 



We investigated properties of DWs in a tetragonal chiral p-wave superconductor. Various cases of the DW structure 
and orientation with respect to in-plane crystallographic axes were analyzed. The magnetic-induction jump on the 
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DW changes from case to case, in contrast to the Landau-Lifshitz theory of ferromagnetism, where this jump has the 
universal value proportional to the local spontaneous magnetization. This conclusion stresses again that the latter is 
not defined for orbital ferromagnetism. The magnetic-induction jump is responsible for stray magnetic fields, which 
have not yet detected outside Sr2Ru04 samples [6]. A quantitative evaluation of stray fields cannot be done without 
a more detailed information (from the microscopical theory or from the experiment) of the parameters of the GL 
theory. 
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